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ABSTRACT 

A functional integral formalism is used to derive the extension of a stiff chain subject to an 
external force. The force versus extension curves are calculated using a meanfield approach in 
which the hard constraint u 2 (s) = 1 is replaced by a global constraint (u 2 (s)) = 1 where u(s) 
is the tangent vector describing the chain and s is the arc length. The theory quantitatively 
reproduces the experimental results for DNA that is subject to a constant force. 

We also treat the problems of a semiflexible chain in a nematic field. In the limit of weak 
nematic field strength our treatment reproduces the exact results for chain expansion parallel 
to the director. When the strength of nematic field is large, a situation in which there are two 
equivalent minima in the free energy, the intrinsically meanfield approach yields incorrect 
results for the dependence of the persistence length on the nematic field. 
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I. Introduction 

Inspired by the elegant experiments by Smith et al. on the response of DNA to 
stetching by a constant force few theoretical papers [2-5] have considered the effect of external 
field on semiflexible chains. Since DNA is stiff it can be described, at least approximately, 
as a worm like chain. The eariest theoretical papers dealing with this problem were initiated 
by Fixman and Kovac JJ]. These authors used a modified version of the Gaussian model 
for stiff chains and provided expressions for the stretching as a function of applied force. 
Their treatment is only valid when the applied force is small and significant deviation from 
these predictions are observed at sufficiently large values of the external force. Marko and 
Siggia [0] have calculated the extension as a function of force for worm like chains and found 
that their results fit the experimental data very well. Some aspects of this theory have 
also been considered by Odijk || who also discusses the competition between entropically 
dominated effects and elasticity effects. 

In this paper we revisit the problem of a semiflexible chain subject to tension. A func- 
tional integral formalism together with a meanfield treatment allows us to set up a general 
way of tackling problems involving semiflexible chains. When our theory is applied to the 
case of stretching of DNA by a constant force, we obtain excellent agreement with experi- 
mental results for the force versus extension. 

E. Elastic Response of a Semiflexible Chain 

(a) Meanfield Theory 

The simplest model for a semiflexible chain (SC) is obtained by taking into account the 
persistence in the tangential direction. A SC can be described by a space curve r(s) of total 
contour length L with s being the arc length. The unit tangent vector is u(s) = dr(s)/ds 
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with the constraint that u 2 (s) = 1 for all s. Since the molecule is stiff it costs energy to bend 
and the bending energy per segment length is proportional to ) 2 - The difficulty with 
calculations involving this formulation of SC is that one has to invoke the constraint that 
u 2 (s) = 1. This makes even the free chain problem non-linear. There has been extensive 
treatment of models of SC in the literature [6-15]. 

Now consider applying an external field that stretches the chain. The effective free energy 
can be written as 

F = -k B T In Z (la) 

where 

Z = J V[u(s)} 5(u 2 (s) - 1) e- n/k * T (lb) 

with 

^4 /( ^ )WfW . uWd , (lc) 

In Eq.(lc) l p is the persistence length of the semiflexible chain (which in the experiments 
of Smith at al. |l] on DNA is estimated to be 53 nm). For generality we have assumed 
that the external field f(s) depends on the arc length s. The case of uniform f(s) = f 
is appropriate for the experiments of Smith et al. and is the one treated in the previous 
theoretical papers. In particular Marko and Siggia || have used an eigenfunction expansion 
to obtain the extension versus / by exploiting the analogy between this problem and the 
quantum problem of a dipolar rotor in an external electric field. 

Here we use a meanfield approach |6|, Q that effectively replaces the local constraint 
u 2 (s) = 1 by a global one (u 2 (s)) = 1. Such a theory has been shown to give an exact 
expression for the mean end-to-end distance of SC when f(s) = |J. The basic idea is 
to enforce the constraint u 2 (s) = 1 using an auxiliary field variable A(s) and evaluate the 
resulting integrals over A(s) by the stationary phase approximation. The free energy (cf. Eq. 
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(1)) can be rewritten as 

/ioo r r rL 

V[X(s)] / V[u(s)] *[u(s),f(s)] exp / X(s)ds 
-ioo J l J0 

/ioo 
\(s)ds exp{— f (s)]} + const. 
-ioo 



+ const. 



where 



and 



with 



*[u(s),f(s)]=expf-i ( L f L u(s) Q(s,s') u(s') + /" u(s) ■ f(s)<L« 

L Jo JO JO 
^[A(a), f (a)] = -In / V[u(s)\ *[u(a), f (a)] - jf* A( S )d S 



(2a) 



(26) 



(2c) 



(9 \2 



Q(s,s') = [-ip(^) +2A(s)]o(s / -s). (2d) 

The constants in Eq. (2a) comes from with normalizations associated with X>[A(s)] and 
2?[u(s)]. These constants will be ommitted from now on. 

The function jF[A(s), f (s)] is a generating functional which can be used to calculate 
various correlation functions. For example the connected correlation function 



(u(s).u(s')) c = (u(s).u(s'))-(u(s))-(u(s>)) 

d dT 
= ~df{s) ' df(s') 

= 3Q-\s,s'). 



(3a) 

(36) 
(3c) 



Similarly AR 2 = (R 2 ) — (R) 2 can be written as 



AR 2 = f L ds[ L ds' 9 QT 



(4) 



lo Jo df(s) df(s')' 

By performing the functional integrations with respect to u(s) the free energy jF[A(s), f(s)] 
can be written as 



.F[A(s),u(s)] = | tr \nQ-±£ds£ds' f(s) Q _1 (s,s') f(s') - J\(s)ds. (5) 



The variable A(s), that enforces the constraint u 2 (s) = 1, is evaluated by stationary phase 
approach || ||. The stationary phase condition is obtained by requiring that be an 
extremum. This leads to the equation 



+ f L ds [ L ds' f(s') Q-\s, s') Q-\s, s") f(s") = 1. (6) 

.QJs,s Jo Jo 

The value of A(s) that makes J-"[X, f(s)] stationary depends on the precise form of f (s). We 
now specialize to the condition that f(s) = f = const, independent of s. This is the situation 
that has been treated in the literature and is assumed to be relevant to the experiments of 
Smith et al. []IJ. For constant f the stationarity condition (Eq. (6)) gives a uniform value 
for A. More precisely Eq. (6) reduces to 
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In order to evaluate the first term in Eq. (7) in a transparant manner we use the boundary 
condition u(0) = u(L) and 9u ^ = dU g^ ■ These conditions were implicitly assumed in Eq. 
(2b). The easiest way to compute the first term in Eq. (7) is in terms of an eigenfunction 
expansion. Let {|s}} denote the eigenstate with s the curviliniear space label and let {\n)} 
be the states that are Fourier conjugate to {\s)} in such a way 

d\n) /2nn\. , , . 

The states |n) are eigenstate of " momentum" such that (n\s) = ^exp( t2 ^ ns ). In terms of 
these eigenstates the first term in Eq. (7) becomes 



n)(n\s) 



1 



2\ds'J ^ A 



>3C 



3 



■2 [ ^ _|_ /p (27m) 2 I" ( 9 ) 



n=— 00 



2 L 



If we use the identity 



^2, cos nx 7r cosh a(ir — x) 1 . . 

^ n 2 + a 2 2a sinh an 2a 2 



the stationarity condition for A becomes 



^ coth(iOL) + i = 1 (106) 



with Q = J^. For large L Eq. (10b) simplifies to 



1 (10c) 

As / — > this stationarity condition coincides with the one derived previously. 



(b) Correlation Function and Mean Square Internal Distance 

In terms of the stationarity solution of A, the various correlation functions can be com- 
puted. For example 

(u(s) -u(s')) = 3Q-\s,s') 4 



f 2 
4A 2 



_ 3 p" cosh[(L-2|^- 3 |)n/2] / 2 

" 4 V/pA sinh(fiL/2) 4A 2 ' 1 J 

For s = s' the above equation reduces exactly to the stationarity condition for A (cf Eq. 
(10b)) leading to (u 2 (s)) = I. Thus our approach ensures that the constraint u 2 (s) = 1 is 
satisfied globally. For L — > oo Eq. (11) becomes 



(u(s)-u(s>)) = ^ 1 ^e^'^ + L (12) 

The mean squared internal distance of the semiflexible chain under tension can be obtained 

as 

<|r( S ')-r( S )| 2 >= [ S 'd Sl [ S 'ds 2 (u( Sl ).u(s 2 )). (13) 

J S J S 

Substituting Eq. (12) in Eq. (13) yields 

<|r( S ')-r( S )| 2 >= r'r'(u( Sl )-u( S2 ))d Sl d S2 

J s J s 
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2 r 2 



\s — s\ 



'I _ e -nk'-sh 



-n 1 1 n 2 

The above equation naturally suggests a length scale, If, below which entropy factors domi- 
nates and above which the mechanical energy associated with the orienting field dominates. 
For small values of the contour length I = \s' — s\ the field dependent term becomes negligible 
whereas for long I the third term in Eq. (14) dominates. A scale If is obtained by balancing 
the first and last terms of Eq. (14) and is given by 

l f = —p-- (15) 

When flf <C 1 then the stationary phase condition can be easily solved and one gets 
A(/) ~ l/l p and consequently Eq. (15) becomes If ~ 1/ f 2 l p . Under these conditions one can 
think of the semiflexible chain to break up into a sequence of "blobs" of effective segment 
length given by 

Z f ~f-\ (16) 

The above result coincides with the blob size predicted by flexible chain under tension. The 
above blob length is called tensile screening length [plj. For Eq. (16) to be valid it is nec- 



essary that within the length £y there should be a large number of segments each of length 
/ p . Since the effect / within £/ is negligible the chain on this length behaves effectively as 
a Gaussian chain containing lf/l p units. It is for this reason in the law / limit coincides 
with the Pincus blob length. In the limit of Iff — > oo, Eq. (15) gives If ~ 1/f and hence 
the chain conformation is dictated by coupling to the mechanical energy. In this limit we 
expect the chain to be aligned with the external field with almost complete suppression of 
fluctuations. The effect of / is felt at all length scales, i.e., If is effectively zero. 



(c) Average Chain Extension 
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The quantity of experimental interest is the average extension z of the chain parallel to 
the external field. The average elongation z is computed using 

2 = <f u(s) 'M>- (17) 

The statistical average in Eq. (17) can be conveniently expressed in terms of the free energy 
functional T. When f is uniform we get 



f SJ 7 , r L 



Ifl 6f J Jo 



5 ^2 



f L "ip 77 +2A(/) 



5s 



i 



ds 



IL (18) 



2A(/) 

In Eq. (18) we have shown the argument of A to emphasize the dependence of A on /. The 
above equation and the associated stationarity condition (cf. Eq. (10c)) determines the 
average chain extension. 

The chain extension z can be easily calculated for the case of l p f — > and for l p f — > oo. 
When lpf — > we get 

* = /f (19) 
where = 2l L = 2(|/ P )L is the size of the corresponding ideal chain with l = |/ p the 
meanfield persistence length M . The above result is also obtained for a Gaussian chain under 
tension. The leading order correction to Eq. (19) can be obtained for l p f <C 1 by expanding 
Eq. (18) in power of /. The stationarity condition for A(/) up to 0(f 2 ) becomes 



P 
2AB 



Kf) ~ A + ^ (20) 



where Ao is the stationary phase condition for / = [BJ. Thus the average elongation becomes 

\ « Ml - Ifll). (21) 

In the small / limit the entropy considerations dominate the effects due to the orienting 
field. 



The dependence of z on / is quite different in the opposite limit. If fl p ~ 1 then we 
expect the external field to be relevant at all length scales. This would suppress the chain 
fluctuations on scale greater than l p . When l p f — > oo the stationary phase condition for A 
has a solution A -1 ~ resulting in A = \f (Cf. Eq. (10c)). In this limit z/L — > 1 and the 
chain assumes a rod conformation. In the limit of l p f > 1 we can approximately solve Eq. 
(10c) to get 

+ (22) 
The chain extension in the limit of large / becomes 




This result has already been noted in the literature ||. The / _1 ^ 2 behaviour for semiflexible 
chains is in contrast to the Gaussian case. Apart from a numercal factor, the result in Eq. 
(23), valid for l p f ^> 1, coincides with that discussed recently by Odijk || for semiflexible 
chains near the rod limit. The prefactor in front of 1/ Jl p f obtained by Odijk is ~ which is 
slightly smaller than our result. The reason may lie in the fact that our treatment utilizes a 
stationary phase method to enforce the hard constraint u 2 (s) = 1 globally. 

In the limit of l p f ^> 1 the chain fluctuations are relatively small and can be expanded 
in terms of Of = (0 X , y ) where Of is the angle between the tangent vector and the external 
field. The average extension is related to (0j(s)) as 

^1-\{0%s)). (24) 

For the uniform field we expect {0"j(s)) to be independent of s hence {0"f(s)) = (Oj). Com- 
paring Eq. (23) we get an estimation of the mean fluctuations 

(«?> - 1/5 < 25 > 

which also coincides with the result of Odijk || apart from the numerical factor. 
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In addition to yielding the results reported in the literature our theory offers a simple 
estimate of the extension z/L as a function of / for arbitrary values of /. This is achieved 
by similtaneously solving Eq. (10c) and Eq. (18). As a test of the utility of the stationary 
phase approach we compare our theory with the experiments of Smith et al. ]I[ on DNA. 
Marko and Siggia have already shown that their numerically exact force extension for 
the semiflexible chain yields excellent agreement with experiments. In Fig. (1) we plot z 
versus / using the parameters appropriate for DNA. For comparison a few points from the 
experiments, as presented in Fig. (1) of Ref. ||, are also shown. The meanfield approach 
reproduces the data quantitatively. 

IE. Conclusions 

In this article we have considered the problem of semiflexible chain subject to an external 
field using a functional integral formalism. The crux of our method hinges on replacing the 
local constraint that u 2 (s) = 1 by a global constraint (u 2 (s)) = 1 where u(s) is the tangent 
vector. The intrinsically meanfield approach has been shown to be successful in producing the 
configurational properties of semiflexible chains in the absence of the external field P, f7L H|. 
Here we have shown that the stationary phase method of enforcing the global constraint 
yields excellent results even when the chain is subject to tension. The very good agreement 
between the theory and the experiments on DNA subject to force is a confirmation of this 
assertion. 

The approach we have described here is sytematic but is not without limitations. These 
limitations become evident by considering the behaviour of a semiflexble chain in a nematic 
environment, a problem that has received considerable attention in the literature [17-22]. A 
meanfield Hamiltonian of a stiff test chain in a matrix of other chains in a nematic state can 
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be written as [19a] 



U = \ [ L [ L dsds'u(s)Q(s, s')u(s') - g T u 2 Js)ds (26) 
Jo Jo Jo 



where u(s) = (u±(s),u z (s)) and g (with the dimension of length) is the strength of the 
nematic potential. By following the procedure outlined in the previous section the stationary 
equation (one that extremizes the free energy with Eq. (26) as the Hamiltonian) becomes 



1_ W + i R(A^)- (27) 
The above equation ensures that 

1 = (ul(s)) + (u z (s)) = (u 2 (s)) 

instead of u 2 (s) = 1 for all s. The mean extention of the chain parallel and perpendicular 
to the director axis can be calculated as 

(R 2 ) = [ L f L (u ± (s) ■ u ± (s'))dsds' (28a) 
Jo Jo 

(R 2 Z ) = f L f L (u z (s) ■ u z (s'))dsds'. (286) 
Jo Jo 

Consider the weak nematic limit, i.e., g — > 0. For small g the stationary condition can be 
solved and the results can be used to get (R±) and (R 2 Z ). These lead to 

27^TT « 1 " Wo (29a) 



j(2Z L) 



\gk (296) 



where Iq = |/ p . The persistence length along the nematic field is increased by a factor of 
(1 + IgZo) which compares well with the exact results (in the limit of L — >• oo) of Warner et 
al. [0. 

A similar analysis for g — ► oo suggests that chain fluctuations perpendicular to the di- 
rector are totally suppressed whereas the persistence length along the director is predicted 
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to increase by a factor of 2. This result is in contrast with the analysis of several au- 
thors [17,19a,21] who have shown that the effective persistence length grows exponentially 
as (g/p) 1//2 . The nematic potential, —gu 2 z , has two deep minima at large g at u z = ±1. Thus 
in the large g limit the chain configuration is dominated by "tunneling " between the two 
minima by instanton [Q. Mathematically the partition function is dominated by instanton 
contributions which our meanfield theory fails to capture. It is clear that if there is symmetry 
breaking in the problem then the replacement of the hard constraint by a global one can 
lead to incorrect results. 
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Figure Caption 

Fig.(l): The solid line is the force-extension curve obtained by simultaneously solving Eq. 
(10c) and Eq. (18). The parameters have been chosen from the fit of force-extension curve 
for DNA reported by Smith et al. |IJ. The value of L = 32.9/im and Iq = 53nm. The square 
represents the experimental results of Smith et al. as reported in Ref. [fll Q . 
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